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We consider a braneworld inflation model driven by the dynamics of a scalar fleld living in the 
5-dimensional bulk, the so-called "bulk inflaton model", and investigate the geometry in the bulk 
and large scale cosmological perturbations on the brane. The bulk gravitational effects on the 
brane are described by a projection of the 5-dimensional Weyl tensor, which we denote by Efj,^. 
Focusing on a tachionic potential model, we take a perturbative approach in the anti-de Sitter 
(AdSs) background with a single de Sitter brane. We first formulate the evolution equations for 
E^u in the bulk. Next, applying them to the case of a spatially homogeneous brane, we obtain two 
different integral expressions for -E^^. One of them reduces to the expression obtained previously 
when evaluated on the brane. The other is a new expression that may be useful for analyzing the 
bulk geometry. Then we consider superhorizon scale cosmological perturbations and evaluate the 
bulk effects onto the brane. In the limit H^li^ <C 1, where H is the Hubble parameter on the brane 
and I is the bulk curvature radius, we find that the effective theory on the brane is identical to the 
4-dimensional Einstein-scalar theory with a simple rescaling of the potential even under the presence 
of inhomogeneities. In particular, it is found that the anticipated non-trivial bulk effect due to the 
spatially anisotropic part of _Ep^ may appear only at 0{H^£'^). 

PACS numbers: PACS: 04.50.+h; 98.80.Cq 



Recent progress in particle physics has lead to the idea that our universe is a (mem)brane in a higher dimensional 
spacetime and all the interactions but gravity are confined on the brane, the so-called braneworld scenario [1]. In 
particular, the scenario proposed by Randall and Sundrum [2], in which a single positive tension brane is embedded 
in the 5-dimensional anti-de Sitter space (AdSs) with Z2 symmetry (the RS2 scenario), has attracted much attention 
as an alternative to the standard 4-dimensional cosmology [3] because of its attractive feature that the standard 

4- dimensional gravity is recovered on the brane in the low energy limit [2, 4]. The inflationary cosmology is also 
widely accepted because of its success in explaining cosmological observations [5] . It is therefore natural to consider 
the inflationary cosmology from the point of view of the braneworld scenario [6]. Brane inflation in the quantum 
cosmological context is also of great interest [7, 8]. 

In this paper, we consider a model of braneworld inflation driven by the dynamics of a scalar field in the 

5- dimensional bulk [9-13] in the context of the RS2 scenario, and investigate the cosmological perturbations on 
superhorizon scales, i.e., in the long wavelength limit. The existence of such a scalar field is supported from the 
unified theoretical point of view, as the reduction of a higher dimensional theory to 5 dimensions will give rise to 
scalar-tensor type effective theories. For a relevant form of the potential in the bulk, the scalar field dynamics can 
give rise to inflation on the brane [9] . A theoretically interesting aspect of this bulk inflaton model is that the bulk 
is not inflating at all; inflation and the subsequent hot Friedmann stage of the universe can be realized solely by the 
geometrical dynamics of the bulk. Furthermore, in the low energy limit H£ <C 1 where H is the Hubble parameter 
of the brane and £ is the curvature radius of the bulk, as far as the spatially homogeneous dynamics is concerned, 
the dynamics on the brane has been found to be indistinguishable form the conventional slow-roll infiation [12]. It 
was also shown that the quantum fluctuations projected on the brane [10] as well as reheating after infiation [13] also 
mimic the 4-dimensional standard infiation, as long as we consider the low energy limit. 

Thus the bulk infiaton model may be a viable alternative to the conventional 4-dimensional inflationary scenario. 
It is then of importance to clarify if the cosmological perturbations produced in the bulk inflaton model also have the 
same desirable features as those in the 4-dimensional theory in the low energy limit, and if there exists any signature 
specific to the braneworld scenario that can be observationally tested. Very recently, cosmological perturbations in a 
bulk inflaton model with an exponential potential and dilatonic coupling to the brane tension have been investigated 
by Koyama and Takahashi [14]. In this paper, we consider a tachionic bulk potential with the potential maximum 
at (j) — 0, with no coupling to the brane tension, and investigate the bulk geometry and cosmological perturbations 



I. INTRODUCTION 
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in this model. We take the geometrical approach to the effective gravitational equations on the brane developed by 
Shiromizu, Maeda and Sasaki [15] and by Maeda and Wands [16]. In this formalism, the bulk gravitational effects on 
the brane are described by the projected 5-dimensional Weyl tensor, 

where is the unit vector normal to the brane (Latin indices run over 5 dimensions, while Greek indices over 4 
dimensions, with the choice of Gaussian normal coordinates such that n^da = 85). 

This paper is organized as follows. In Sec. II, we briefly review the bulk inflaton model, based on the geometrical 
formalism [15, 16]. In Sec. Ill, we derive the evolution equations for E^i, in the bulk when there exists a non-trivial 
bulk energy momentum tensor, and solve them in the case of a spatially homogeneous brane, under the assumption 
that the amplitude of the scalar field (f> is small. In Sec. IV, we consider the cosmological perturbations on superhorizon 
scales on the brane by solving the perturbation equations for (f> and E^^ in the bulk. Since the tensor perturbations in 
this model are identical to the vacuum AdS^ bulk model at leading order, we focus on the scalar-type perturbations. 
We find the standard 4-dimensional result is recovered in the low energy limit iJ^f^ <^ 1, and possible signatures of 
the braneworld appear only at 0{H^£^). In Sec. V , we summarize our results and discuss the implications. Some 
formulas used in the text are summarized in Appendices A and B. Properties of the Green function for i?^^ in the 
bulk are analyzed in Appendix C. 



II. REVIEW OF THE BULK INFLATON MODEL 

First, we review the bulk inflaton model [9]. We choose the Gaussian normal coordinates, 

ds"^ = {nam + qab)dx°'dx^ = dr^ + q^^dx^dx" . (2.1) 

We assume that the brane is located at r = ro . By extremizing the action of the 5-dimensional Einstein-scalar system 
with a brane, the gravitational equations in the bulk take the form 

^^^Gab + A59ab = 1^1 ( Tab [</>] + ^a6<5(r " ro)) , (2.2) 

where and A5 are the 5-dimensional gravitational and cosmological constants, respectively. Tab [4>] is the energy- 
momentum tensor of the bulk scalar field for which we assume the form. 

Tab M = dact>db<l^ - gab(l9''''d,(t>dd(l) + y [</)]) , (2.3) 
and Sab is the energy-momentum tensor on the brane for which we assume the vacuum form. 

Sab = -CTQab, (2.4) 

where a is the brane tension. Following the spirit of the RS2 scenario, we assume the tuning between A5 and a as 



A. = -^ = -ii, (2^5) 

where £0 is the AdSs curvature radius of the RS2 braneworld, and assume the Z2 symmetry with respect to the brane. 
Furthermore, we assume that the potential takes a tachionic form 

y(<^) = + ; 1^0 >0, m2<0, (2.6) 

in the vicinity of (/) ~ 0. We tacitly assume that there exists a minimum of the potential somewhere at (p = 4>min 7^ 
at which V{(t>ynin) = 0, where the RS2 flat brane is recovered. 
The field equation for in the bulk is 

V"V„^-^1/[0]=O, (2.7) 
where Va is the 5-dimensional covariant derivative. Here, for simplicity, we impose the Neumann boundary condition 

dr(t)\r=ro = , (2.8) 

on the scalar field. This implies that the scalar field does not couple to the brane tension. 
From Israel's junction condition on the brane. 



-uli^S^^ - ^Qi^-yS^ = -^Klcrq^cu , (2.9) 
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and the Zi symmetry, we obtain the 4-dimensional effective gravitational equations on the brane [9, 15] 



(4) 



where 



(2.10) 

(2.11) 
(2.12) 



'^^^Cabcd is fivc-dimcnsional Weyl tensor. 

In general, the projected Weyl tensor term E^,j cannot be determined without solving the bulk dynamics. However, 
for a spatially homogeneous and isotropic brane, i^^^ on the brane may be evaluated without solving the bulk. By 
using the 4-dimensional contracted Bianchi identities, one finds [9] 



(2.13) 



where a{t) is the cosmic scale factor on the brane and t the cosmic proper time. The other components of E^i, are 
determined by the isotropy of the brane and the traceless nature of E^i,. 

For the potential of the form (2.6), we then focus on the region \rn?\(j)'^ <C Vb and consider the perturbation with 
respect to the amplitude of <f). In the zeroth order, the background metric is determined as [7] 



ds^ = b^{z){dz^ - dt^ + a{tf5ijdx'dx^), 



(2.14) 



biz) 



HI 



ait) 



sinh[iJ(|2| +^o)] ' H 
Here dz = dr/b{r) and the brane is located at 2; = 0, H£ = sinh[if^o]) and 



A5 + 



is the effective AdSs radius. Note that i > £0. The braneworld inflates with the Hubble rate H given by 

Vo 



(2.15) 



(2.16) 



In the first order in (j), the evolution of the scalar field on the background metric is determined. Expressing (f> as 
the sum over all possible modes, (j) = Un{z)tpn{t) , the late time behavior of </> on the brane is dominated by the 
bound-state mode for which the effective 4-dimensional mass-squared for (j){t) is given by [9] 



m 
~2 



(2.17) 



for H^r <C 1 and \nr\/H^ <C 1. In this case, one has 



m 

T 



(2.18) 



on the brane. 

The effective equations (2.10) at second order is then derived as follows. Using Eq. (2.18), Eu is evaluated as 



E, 



J 



+ 2H4>) = - 



^5 12 



+ 



c 



(2.19) 



where C is an integration constant which depends on the initial condition. The term C/a* is called dark radiation 
term. In the present case, since we are interested in the late time behavior at the inflationary stage, the dark radiation 
term can be safely neglected. Then, the effective Friedmann equation becomes 



3(-) =kIp^''^ -E„ = Klpes., 



where 



K4 



4 ' 



(2.20) 



(2.21) 
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and 

Peff = 4 Q-^' + ^^(0)) • (2.22) 
Therefore, by introducing the effective 4- dimensional field ^p as 

^=A</'. (2-23) 
the system is equivalent to the 4-dimensional Einstein-scalar system with the potential, 

Viip) = . (2.24) 

Thus the effective field (p behaves as a conventional inflaton on the brane, and the conventional slow-roll inflation is 
realized on the brane in the low energy approximation H^£^ <C f [9, 12]. 



III. EVOLUTION EQUATIONS IN THE BULK 

When the bulk scalar field is spatially inhomogeneous, it is necessary to solve the bulk geometry to evaluate E^^ 
on the brane. In this section, we derive the evolution equations for E^,^. We then apply the resulting equations to 
the spatially homogeneous and isotropic case to see how the previous result reviewed in Sec. II is recovered from the 
bulk point of view. 



A. General equations 

First, we recapitulate the definitions of the two projected Weyl tensors in the bulk, 

Eab = Cced/n^'n'^ ql ql , Babe = ql ql CdecfU* . (3.1) 

By definition, Eab is symmetric with respect to the indices (a, b) whereas Babe is anti-symmetric with respect to (a, b). 
Furthermore Eab is traceless E^a = 0. 

We start from the 5-dimensional Bianchi identities 



(5) 

V[o Rbc]de = 0, 



(3.2) 



(5) 



where Rabcd is the 5-dimensional Riemann tensor. We can derive the equations for E/j^^ and Babe from Eq. (3.2) by 
using the 5-dimensional Einstein equations (2.2). We find 



where 



with 



(4) 3 1 

£nEab = -C^-Bc(ah) + Cacbd + 4i^(^„£^(,)c " -^KEab — ■^qab^'^ Eed 

+2KaeK'"'Kbd - ^-KabK^'^Ked + 4 Pab [Tab] 



Pab ^Pab 



■^Tden'^rfqabK + ^TKab - \Kf^ql)Tde + Ik'^'Tm 

+ l£nFab - '^D-(Taen^q(^q(^a]b)) - D(^a[qt)Tden''] , 



(3.3) 



(3.4) 



and 



where 



Fab = Tf.qiql - {Tfgnfn^ - -T)qab - ^Tf^qf^qab, 



Kab = Kab - -Kqab , 



£nBabe = -2-D[(j£^(,]c "I" K^Babg — 2Bcg[aK^. + Kg Qabe[Tab] 



_ 2 ^ ^ \ 2 

Qabc[Tab\ = --^Tden''[q^a^b]e- Ky^qb]cJ + 



Tden''q[aqb]c 



qb]c{Tden'^n'' 



:D,, 



it^llTde 



(3.5) 
(3.6) 



(3.7) 



They are first-order partial differential equations for Eab and Babe- 
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B. Second order equations for Ef^v 



We assume the tachionic potential form (2.6), and take the perturbation approach with respect to the amplitude 
of (f). 

First we expand the energy momentum tensor of the scalar field as 



(3.8) 



So far, we have not specified the coordinate system. Now, for definiteness, we choose the conformal (Gaussian normal) 
coordinate (2.14). The extrinsic curvature of each z = constant hypersurface is expanded as 



26(z) 
4? = 



(3.9) 



corresponding to the amplitude of the scalar field. 

Substituting the above expansions of Tab and Kfj,^, into Eqs. (3.3) and (3.6), eliminating B^„a in favor of E^^^, and 
keeping terms up to 0{(j)^), we obtain the desired second order equation for as 



where 



4, = h''{z)E^,, 

C^b{z)^-l-J-+n,-4H\ 

oz b[z) oz 

SAT!t^]^b\z)[D-Q^,,^)[T^^]+£^[p,^[T^^^^ 



(3.10) 



(2) 



(3.11) 



where Tia^-* is the second order part of t!)^ given by Eq. (2.11). The operator C is the same as the one derived 
in [17]. What is new is that the source term due to the presence of the bulk scalar field is taken into account. The 
boundary condition on the brane is 

d_ 



dz i^'^") b ~ '^^ ^^'^ab] I 



(3.12) 



which is obtained from Eq. (3.3) and Israel's junction conditions (2.9). Thus, E^^^ has a jump at the brane in general, 
whereas (j) is smooth as given by Eq. (2.8). 

In addition to the above, we have a set of constraint equations that come from the 4-dimensional contracted Bianchi 
identities on each z =constant hypersurface. They are 



d^e^, = kIyi,[tI 



(2)1 



(3.13) 



If we regard E^^ as the 'electric field', these constraints imply that S,^[T^^^] plays the role of the 'charge density' in 
the bulk. 

Although the bulk equation (3.10) is perfectly legitimate, we find it is more convenient to deal with equations that 
are reduced from Eq. (3.10) by using the constraint equation (3.13), which we call the reduced equations. 
The components of Eq. (3.10) are explicitly written down as 



Kz)dJ^^d,)Eu + ^H^Ett - d^Eu - SHdtEtt + ^^^Eu - ^O'^Ekt 
\b{z) I 



'2U 

2HdiEtt 7rd Eki 



1^5 ^tt, 



'^5 Sti , 



(3.14) 



(3.15) 



b{z)d, - d^Eij + HdtEij 



+^V^Eij - 2H{diEtj + djEu) + 2a^H^6ijEt, 



^5 ^ij ■ 



(3.16) 
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In the above equations, the components of arc coupled each other, which makes it very difficult to solve them. 

However, this situation can be improved by using the constraint equation (3.13). The components of Eq. (3.13) are 
written down as 

-dtEu + ^d'^Ekt - ^HEtt = 4 b^{z)j:c , (3.17) 
-dtEu + -^d^Eki - 2,HEu = 4 b''{z)Y.i . (3.18) 
Using these, we may eliminate d'^Ekt and d'^Eki from Eqs. (3.14) and (3.15), respectively, to obtain 

b{z)d, (^a,)4t - UH^Ett - d^Ett - IHdtEtt + ^V^^tt = ^ Su ; 

Stt = Stt + AHb\z)^t . (3.19) 
b{z)d,(^d,)Eu - 2H^Eu - d^Eu - SHdtEu + ^V^E^^ = 4 Su ; 

Su = Su + 2Hb\z)i:i + ^d,Ett, (3.20) 

and Eq. (3.16) is rewritten as 

b{z)d, (^-^d.^Ei, - d^Eij + HdtEij + ^V^iy = kI Si, ; 

Sij = Sij + \ (2H{diEtj + djEu) - 2a^H^5ijEu) • (3.21) 

Thus, one can determine the components of E^i, step by step; solving first Eq. (3.19) for Eu, next Eq. (3.20) for Eu, 
and finally Eq. (3.21) for . 



C. The spatially homogeneous and isotropic background 



Here, we solve the evolution equation (3.19) in the spatially homogeneous and isotropic case. In this case, we have 
already seen that E^i, on the brane can be obtained without solving the bulk, as given by Eq. (2.19). The purpose 
of this subsection is to determine E^^, in the bulk and to take the limit of it to the brane to examine if E^^ thus 
obtained gives the correct answer for E^,^ on the brane, as a check of our derivation given above. 

Equation (3.19) in the spatially homogeneous case is 



6(^)9.(^5.) - - 7Hdt - 12H' 



Et, 



Kr St 



5 '^tt ■ 



To solve this, we introduce the retarded Green function that satisfies 

b{z)d, (^5^) - - 7Hdt - 12if2j G{t, z; t', z') = -6{t - t')6{z - z')b{z). 

The Green function satisfies the reciprocity condition, 

G{t,z;t',z')=G{-t',z';^t,z), 
and the Neumann- type boundary condition on the brane, 

(9^G(t,0+;t',z') = 0. 

An explicit construction of the retarded Green function is given in Appendix C. 
Using the Green function, the formal solution of Eu is expressed as 



Eu{t,z) 



■ / dt'—-G{t,z;t',z')d,^Euit',z') 
Jo b{z') 



r°° dz' 

Jo b{z') 



z'=0 
I J 



dt'G{t, z; t', z') - G{t, z; t', z')dt' + 7HG{t, z; t', z') Eu{t', z') 



fOO pOO 1 

-4 dt' dz'—-G{t,z;t',z')Stt{t',z'), 
Jo Jo b(z') 



t'=0 



(3.22) 

(3.23) 

(3.24) 
(3.25) 



(3.26) 
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where we have assumed a regular initial data at f = so that there is no contribution from z = oo. In the above, the 
first term is the contribution from the brane-boundary of the extra-dimension, the second term is that from the initial 
surface, and the third term is that from the source in the bulk generated from the scalar field. In Appendix C, we 
show that the Green function G{t, z; t', z') at late times behaves as oc a{t)~'^. Therefore, the second term behaves as 
radiation, i.e., it describes the dark radiation term. Thus for a de Sitter brane, or for any observer on the z =constant 
hypersurface, we can neglect it after a sufficient lapse of time. 

The third term can be evaluated as follows. Using the field equation (2.7), we find the relations among several 
terms in the source term as 



£)«Q„« + -St = 0, 
a 

In particular, the latter equality implies the existence of an integrability condition. Namely, if we define 

fit, z) = a\t)b\z)Ptt , g{t, z) = -a\t)b\z)^t , 
Eq. (3.28) implies the existence of a function M{t, z) such that 

fit, z) = d,M{t, z) , g{t, z) = dtM{t, z). 
Introducing a rescaled function of M{t, z) as 

M{t,z) = a-'^{t)M{t,z), 

the source function Su is expressed as 



Stt{t,z) = [&(z)a,(^5,) - - 7Hdt - UH^'\M{t,z). 



(3.27) 

(3.28) 

(3.29) 
(3.30) 
(3.31) 

(3.32) 



We see that the operator acting on M{t, z) is the same as that acting on the Green function in Eq. (3.23). Therefore, 
by repeating integration by part twice, the third term of Eq. (3.26) can be rewritten in the form. 



/ dt'—-G{t,z;t',z')d,,M{t',z') 
Jo o[z ) 



■M{t,z). 



(3.33) 



z'=0 



Here again, assuming the regularity of the scalar field, we have neglected the contribution from z = oo. The boundary 
condition of M{t, z) on the brane-boundary is obtained from its definition as 



d,M{t,(i+) = Ptt{t,(i+). 



(3.34) 



This is the same as the boundary condition for Eu, Eq. (3.12). Thus the boundary contribution in Eq. (3.33) just 
cancels the first term in Eq. (3.26), and we have 



Eu{t, z) = 4M{t,z). 



(3.35) 



that is, the solution for is given by M , apart from the initial surface contribution that can be neglected at late 
times. 

From Eqs. (3.30) and (3.31), it immediately follows that two different representations of M, hence of Ett, are 
possible; the f-integral representation and the ^-integral representation. 

1. f-integral representation 

The f-integral representation of the solution is 



Taking the limit of it to the brane, and eliminating the mass term using the field equation, we obtain 

^2 /•* 



Eu{t,0+) 



1 



a^{t) 



dt'a''{t'){d^,cl>+^^)^ + a\0)E„{0,0+) 



(3.36) 



(3.37) 



where we have recovered the initial surface contribution. This is the same as Eq.(2.19) obtained from the 
contracted Bianchi identities on the brane. 
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2. z-integral representation 

The z-integral representation of the solution is 



Ettit,z) = 



1 



b{zy 



4 I dz'b\z')Ptt{t,z')+Ett{t,0+) 
Jo 



(3.38) 



This is a new form that has not been obtained before. A nice feature of this representation is that it automatically 
satisfies the boundary condition on the brane. One immediate consequence is that the regularity at 2; = oo implies 
the equality, 



POO 

Ett{t,0+) = -4 / dzb^{z)Pt, 
Jo 



{t,z). 



(3.39) 



Thus, Ett on the brane may be expressed in terms of an integral of over the extra-dimension on the t =constant 
hypersurface with the integrand given by a certain combination of second derivatives of the energy momentum 
tensor. A point to be emphasized is that Eq. (3.39) has no explicit dependence on the initial condition; Ett on 
the brane is solely determined by the field configuration on the t = constant hypersurface. This is similar to the 
local mass conservation law that holds in a spherically symmetric spacetime. Apparently, Eqs. (3.37) and (3.39) 
must coincide with each other up to the term proportional to a~'^{t). This is guaranteed by Eq. (3.28). The 
representation (3.38) may be useful for analyzing the behavior of E^f in the bulk, though we have not explored 
it yet. 

Given the solution for Ett, it is trivial to obtain the other components Eti and Eij. For Eti and the traceless part 
of Eij, the source terms for a spatially homogeneous (j) are easily found to be zero. Hence they are zero except for the 
contributions from the initial data that decay rapidly in time anyway. The trace part of Eij is simply given by Ett 
because of the traceless nature of E^,^. 



IV. COSMOLOGICAL PERTURBATIONS 

To investigate cosmological perturbations on the brane, we take the following approach. First we derive the 
perturbation equations for the scalar field Scj) and the projected Weyl tensor in the bulk. Then we solve them 
with the assumption that the scale of interest is much greater than the Hubble horizon size, and take their projections 
on the brane to evaluate their effects. 

Note that there are two small parameters in our approach. As before, we assume the amplitude of </> to be small. 
In addition to it, we assume S(t> is much smaller than cp, say S(j) — 0{e(j)) with e ^ 1. In what follows, we consider 
the perturbations accurate to 0{4P') in the scalar field amplitude and linear in e. An important consequence is that 
the perturbation of the energy momentum tensor becomes effectively gauge-invariant as will be shown below. Our 
calculations are considerably simplified by this fact. 

We focus on the so-called scalar-type cosmological perturbations, because the tensor-type perturbations (which are 
spatially transverse-traceless on the brane) are identical to those in the vacuum AdSs bulk model discussed in the 
literature [7, 18] to the accuracy of 0{(j)'^). 



A. The perturbation equations in the bulk 



First, we write down the perturbation equations in the bulk. As usual, we expand all the perturbation variables in 
terms of the spatial scalar harmonics Y{x^) that satisfy 



{5'^didj + k'^)Y = Q, 
and the associated vector and tensor harmonics 



Yi = ~diY, Yij 



l.d,d,Y +^-5,jY , 



(4.1) 



(4.2) 



where we have assumed spatially flat slicing of each z =constant hypersurface for simplicity. 
The perturbation of the scalar field equation is simply given by 



^d,b^d,-\dta^dt-%-m%-' 
0'' o'' a'^ 



X{t, z) = Q, 



(4.3) 



with the boundary condition dzx\b = 0, where we have set 5(j) = x{i^z)Y{x^). Note that we do not need to take 
account of the metric perturbation in the bulk because it will be of 0{(f)^). 



To write down the perturbation equations for 6Eui,, we expand 6Eui, in terms of the scalar harmonics as 



5Eu = EY , 
5Eu = aEiYi , 

SEij = a^^EYSij + E2Yij) , 



(4.4) 



and their source terms S^^IST^^^] and T,,^[dT^l^] as 

Su[ST;^1^]=SY, 



Su[5Til^] = aSrYi, 

Sij[ST^^] = a^lsY6ij+S2Y,j) 



E,[ST^f] = EY, 
i:i[6T^l^] = ai:^Yi 

The evolution equations of SE^j^t, are written down as 



(4.5) 
(4.6) 



E = kIS + AHb'^T, , 
El = kI {Si + 2i?6^Si) - 2-HE , 



bd,b-^d, - a-^dta^dt - 2ij2 _ _ 



£2 = 482- 4^-HEi , 



(4.7) 



and the constraint equations are written as 

- (dt +4H^E+ -El = kI b*J: 
k 



■ (dt + 4H^Ei -^{E- 2E2) = 4 fe^Ei 



(4.8) 



Using the expression of Sfj^^[6T^^^'] given in Eq. (3.11), we may further deduce the source term as follows. We 
decompose P^i^lST^^^] as 



PY, 



Pij[ST^^]=a'{lPY6ij+P2Yij), 



(4.9) 



andD«Q„(^,)[JT^^)] as 

D'^Q^ulST^lh ^ QY, 

D"Q,,^,j)[ST^f] = a^QoYdij+Q2Y,,), 
where, from the traceless condition of the whole source term, Qo is expressed as 



go = 3 



Q + a-^dt{a^i:) - 



(4.10) 



(4.11) 



With these expressions, S, Si and 5*2 arc expressed as 

S = b''[b-^d,{b''P) + Q + dtE'\, 

Si = b^[b-^d,{b^Pi) + - 1{^S - adt{a-^Ei)} 



S2=b* 



(4.12) 
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The explicit forms of these source terms in terms of the scalar field perturbation are given in Appendix B. 
The boundary conditions of 6E^^ on the brane are 

d,E = 4P, d,E^ = nlPu d,E2 = KlP2, (4.13) 

by using Eq. (3.12). 

The bulk effects on the brane can be analyzed by solving first the scalar field perturbation equation (4.3), inserting 
it to the source term of Eq. (4.7), and solving it under the boundary conditions (4.13). Although to solve them 
analytically is quite difficult or almost impossible, it is possible to solve them in the long wavelength limit k/aH <C 1. 



B. Scf) in the long wavelength limit 
Let us analyze the behavior of the scalar field perturbation x in the long wavelength limit. We expand it as 

/"OO 

X{t, z) = Xb{t, z) + dMxM{t,z), (4.14) 

J3H/2 

where the first term denotes the contribution from the bound state mode with the effective 4-dimensional mass- 
squared Mgjj = (under the assumption <C H^), and the second term is a superposition of the continuous 

Kaluza-Klein (KK) niodc;s. Since the KK modes have mass-squared larger than 9if^/4, once their wavelengths 

become larger than the Hubble horizon size, they will damp out exponentially as a~^/^. Hence we expect that their 
contributions to be unimportant. At least in the limit H£ <^ 1, there are many pieces of evidence that support this 
expectation. As in the conventional inflationary scenario, the scalar field perturbations will come from the quantum 
fluctuations, and this expectation should be carefully examined in the context of quantum theory. Here, however, we 
simply assume so and focus on the bound state mode. Then we have the relation, 

^d,{b^d,x) = [m^^ -Hpjx^O. (4.15) 

In particular, on the brane, this gives 

d'A = %X = -]-^^^{c?^^x) - • (4.16) 



(4.17) 



C. SE^i, in the long wavelength limit 

By inspecting the explicit forms of the source terms S, Si, S2, S and Si, we flnd their behaviors as 

5, S = 0(1), 5i , El = 0{k), S2 = 0{k^) for ^ 0. 
Then from Eqs. (4.7) we find 

E^0{1), Ei=0{k), E2^0{k^). (4.18) 

Turning to the constraint equations (4.8), we then sec that the long wavelength limit of E and Ei can be obtained 
by simply integrating Eqs. (4.8), i.e., without solving Eqs. (4.7), while we need to solve the 2nd order bulk equation 
to obtain E2. Therefore, let us first solve for E and Ei. Discussion on E2 is deferred to the next subsection. 
At leading order in the long wavelength limit, Eqs. (4.8) reduce to 



Thus E is given by 



where 



-^dt{a''E) = Klb''i: + 0{k^), 
-l^dtia^Ei) = 4 ft^Ei + ^E + 0{k% 

CI oCt 



E{t, z) = -kI ^ y * dt'a\t')^{t', z) + 0{k^) , 



E = 



262 



h' 



i^'x + h') + {<l>'x: + ^'x') - {<l>"x + h") - 2H^X 



(4.19) 
(4.20) 

(4.21) 
(4.22) 
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With thus given E, Ei is given by 



where 



Si = - 



262a 



(4.23) 



(4.24) 



On the brane, the time integrals in Eqs. (4.21) and (4.23) can be exphcitly performed. Using Eq. (4.16) and the 
boundary condition </>' = x' = 0, S on the brane can be rewritten as 



Hence we obtain 



E{t,0+) = -^(l>X 



(4.25) 



(4.26) 



where C is an integration constant. Since the term C/a* decays rapidly, we may neglect it at late times. Then, with 
a similar manipulation, we find 



which gives 



o 0/ 



(4.27) 



(4.28) 



where we have neglected the term oc a ^ as before. Thus we have obtained E and Ei on the brane on superhorizon 
scales. 



D. Bulk anisotropic effect 



Now we turn to solving E2 that describes the spatially anisotropic part of SE^i,. In the long wavelength limit, the 
equation for E2 in Eqs. (4.7) reduces to 



E2=KiS2-4-HEi+0{k^). 



(4.29) 



The El term on the right-hand side of the above equation is given by the time integral in Eq. (4.23). In general, 
it is impossible to perform the integral explicitly. However, if we assume that (f) is dominated by the bound-state 
mode, its time variation is small; i.e., we may adopt the slow-roll approximation. We have already assumed that the 
perturbation x is dominated by the bound-state mode as well. Hence, under the assumption <^ H"^, Ei to the 
lowest order in is found as 



6 aH 



'X + H^x - -^<t>'x 



(4.30) 



Note, however, that this approximation will not be valid at large values of H\z\, since the regularity of the spacetime 

at H\z\ = 00 implies the breakdown of the bound-state mode dominance. Nevertheless, it seems physically reasonable 
to assume that the precise behavior of the source term of E2 at large values of H\z\ will not significantly affect the 
value of E2 on the brane. Under the same approximation, S2 is given by 



S2^-j^(^"x + 2H,Px) . 



Prom Eqs. (4.30) and (4.31), the source term is evaluated as 

o 262 fc2 



KiS2 



AH-Ei 
a 



<t>"x + T^'X 



(4.31) 



(4.32) 
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Given the source term, the formal solution of E2 is given by 

E2{t,z) = - dt' dz'—-G2{t,z-t',z')UlS2-AH-EA{t',z'), 
Jo Jo o[z') V a / 

where G2 is the retarded Green function for E2 that satisfies 

bd^b-^d^ - a-^dta^dt - 2H^]G2{t, z;t', z') = -b{z)S{t - t')S{z - z') , 



(4.33) 



(4.34) 



in the long wavelength limit. In the above, we have neglected the contribution from the initial surface t = since 
is assumed to be very near the top of the potential at t = 0. The contribution from the brane-boundary is absent 
because P2 = on the brane. For JfZ <C 1, and at late times, the retarded Green function behaves as 



G2{t, z; t', z') = {Hl)e{t - t') a-^{t) a^{t') , 

as shown in Appendix C. 

Then, at late times, the ^-integral in Eq. (4.33) can be easily evaluated. From Eq. (4.32) we find 



(4.35) 



f 

Jo 



dz 



1 

W) 



k \ r°° 

]S2-4H-Ei){t,z) oc / dzibcj}" + b'4>')x 
a ' Jo 

"(60'X)' - 



dz 



dzl 



6Y. 



(4.36) 



Now, imder the assumptions that (f) and x ^ire both dominated by the boimd-state mode and m^C'^ <C H^f^ ^ 1, the 
derivatives cf)' and x' s-re both of 0{m?). Hence E2 turns out to be of O(m^). This is a very interesting result. As 
long as we focus on the efi'ects that persist after a sufficiently long lapse of time, signatures of the braneworld in this 
bulk inflaton model do not appear at O(rn^) but only at O(m^). Since we have neglected O(m^) terms when deriving 
the source term (4.32), we are not able to evaluate the O(m^) corrections quantitatively here, but leave it for future 
work. It may be worth mentioning that analyzes of the decay of the scalar field out to the bulk were carried out 
previously and this effect was shown to appear at the same order, i.e., at O(m^) [12, 19, 20]. 



E. Cosmological perturbations on the brane 

We are now ready to discuss cosmological perturbations on the brane. General behavior of superhorizon scale 
cosmological perturbations on the brane has been analyzed with the geometrical approach in [21]. Here we follow 
their analysis. Since the scalar field perturbation induces only scalar- type perturbations on the brane, the perturbed 
metric on the brane 7^1^ = q^vlb^{^+) may be expressed as 

7« = -(i + 2^y), 

7ti = -a{t)BY,, 

7ij = a" {t) {5,j + 2Hl Y5,j + 2Ht ) . (4.37) 
The effective energy-momentum tensor on the brane is 

Tjr^^=io(T^J-^E,,y (4.38) 
where T^xi) is given by Eq. (2.11). The perturbed effective energy-momentum tensor JT^t*^^ is denoted by 

ST^''''^\ = qYi, 

^j.(eff)i^ =p(^T,^Y5i+UYj), (4.39) 

where 5, q, ttl and H represent the density contrast, momentum density fiuctuation, isotropic pressure perturbation 
and anisotropic stress perturbation, respectively. By taking the perturbation of T^^^ and E^^, explicitly, these matter 
perturbation variables are expressed in terms of 5^ and as 



,2k ■ 1 
^3a' 
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5-1 1 



(4.40) 



where it is understood that all the quantities are those evaluated on the brane. 

We have just seen in subsection D that E2 appears only at O(m^). Since m'^/H^ <C 1 is assumed, this implies that 
anisotropic stress is of 0{m^£^) in the low energy approximation H^£^ ^ 1. Hence, to the leading order in the low 
energy expansion, the anisotropic stress perturbation vanishes; 11 = 0. The other matter perturbation variables are 
expressed in terms of x by substituting to the above the expressions of E and Ei, Eqs. (4.26) and (4.28), respectively. 
We thus find 

pS = io{^X+lrn^^x), (4-41) 

q = £o-<j>x, (4.42) 
a 

PTTi, =4(^X- ^m^^x), (4.43) 

in the long wavelength limit. A notable fact is that those expressions are exactly those one would obtain for a 
4-dimensional scalar field with mass-squared M^g = rv? 12, with a simple rescaling ^ if we neglect the 

metric perturbation in them, just as in the case of the homogeneous background. Thus, provided that we can justify 
the neglect of the metric perturbation in these expressions, the description in terms of the effective 4-dimensional 
theory with mass-squared m'^/2 turns out to be valid even for inhomogeneous field configurations, at least for those 
inhomogeneities whose scales are larger than the Hubble horizon scale and when the low energy approximation is 
valid. 

Now, let us justify our neglect of the metric perturbation in the matter variables. The gauge transformation of the 
metric perturbation hjj,i, induced by an infinitesimal coordination transformation = ^ ^'^ is given by 

= h^^ - C^i^ - C^i^ , (4.44) 

where the vertical bar | denotes the covariant differentiation with respect to the 4-mctric 7^,^. Let us assume that a 
certain 'geometrical gauge' is chosen when wc solve the effective 4-dimensional Einstein equations. Here a 'geometrical 
gauge' is a gauge for which the gauge condition is given purely in terms of geometrical quantities. An example is the 
spatially flat slicing TZ = Hl + Ht/3 = or the shear- free (Newton) slicing kB — H'j, = 0. In such a gauge, since the 
metric perturbation is due to the scalar field perturbation we have /i^j^ = 0(e (jy^), where e stands for the amplitude 
of X relative to 4>. Then, for a gauge transformation between two different geometrical gauges, we have = 0{e(f>^). 
This implies that the change of x by such a gauge transformation is x — x = -^<^|^^'^ = 0{e (j)^). Similarly, the change 
of ST^^ is of 0(e (/)*) while ^T^JJ is of 0{e (fp). Hence, to the accuracy of 0{(j)'^) that we are interested in, x as well as 
(5r®^ are gauge-invariant if we restrict our choice of gauge to a geometrical gauge [22] . 

The important point to be kept in mind is that the perturbation equations we derived are valid only in geometrical 
gauges. This is because, if we consider a gauge whose condition involves matter variables, the gauge transformation 
from a geometrical gauge to such a matter-based gauge will give = 0(e), and hence h^i, = 0(e). But our equations 
can produce of O(e0^) only. 

Nevertheless, we may of course introduce various gauge-invariant quantities that may be defined in matter-based 
gauges, as in the standard cosmological perturbation theory [23, 24]. To repeat again, what we should keep in mind is 
that we should evaluate them in a geometrical gauge. As we see below, the most convenient choice is the flat slicing. 
That is, we simply regard our perturbation equations for x as those on the flat slicing. 

Among various gauge- invariant quantities, a convenient one for our discussion is the curvature perturbation on the 
homogeneous density hypersurface, first introduced by Wands et al. [25], 

C-H, + ^ + -J^=n+-J^^. (4.45) 
3 3(p + p) 3{p+P) 

The advantage of using this quantity is that it is known to be conserved in time on supcrhorizon scales when the 
entropy perturbation is negligible, irrespective of the gravitational theory. Since our T^^ is not the real energy- 
momentum tensor in the 4-dimensional Einstein theory, using ( is preferable. Then, choosing the flat slicing TZ = 0, 
from Eq. (4.41) and the fact that p+p = £0^'^, we flnd 

X + M^gCP X ^, M^gCl) X ^, gX .^^g. 

where we have used the slow-roll equation, 3H^ w ~^efi^- "^^^^ same as the one we would obtain for the 

standard 4-dimensional theory. 
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Alternatively, we may choose the curvature perturbation on the comoving hypersurface defined by [24, 26] 



Tic =11- 




(4.47) 



Provided that T°|j may be regarded as the actual 4-dimcnsional energy-momentum tensor, which is in(lcc;d appropri- 
ate [21], TZc is also known to be conversed and equal to C with good accuracy on superhorizon scales. Again, simply 
choosing the flat slicing 7?. = 0, we see that the result is in agreement with Eq. (4.46) within the accuracy of our 
interest. 

Thus, we conclude that all the predictions we obtain in this bulk inflaton model are indistinguishable from the 

case of the 4-dimcnsional inflaton model with mass = ?77.^/2, including the perturbation, in the limit H'^i'^ < 
1. In particular, one would obtain the standard, almost scale-invariant spectrum for the large scale cosmological 
perturbations. In this model, possible signatures of the braneworld appears at most at 0{m^i^), or at second order 
in the slow-roll parameter rn? /H"^. 



V. CONCLUSION 

We investigated superhorizon scale cosmological perturbations in the brane inflation model in which slow-roll 
inflation on the brane is induced by the dynamics of a scalar field (f) living in the bulk. We took the geometrical 

approach in which the projected Weyl tensor E^^ describes the gravitational effect on the brane from the bulk [15]. 

First we derived the evolution equations for E^j^i, in the bulk by assuming a tachionic potential for (p and and focusing 
on the dynamics near the top of the potential with negative mass-squared < 0. We applied them to the case of 
spatially homogeneous and isotropic background. We found two difl'erent integral forms for iJ^y, one with respect to 
time t and the other to the extra-dimensional coordinate z. By taking the brane limit of the f-integral form of E^^i,, 
we recovered the resiilts that had been previously obtained without solving the bulk [9]. Namely, in the low energy 
limit H^£^, by a rescaling of (p to an effective 4-dimensional field, the dynamics on the brane is indistinguishable from 
the standard 4-dimensional infiaton model with the same tachionic potential but multiplied by 1/2. In addition, the 
2-integral form of E^^, which may be useful for analyzing the bulk geometry, is found to give a new expression of E^^, 
on the brane that depends only on the field configuration on the t =constant hypersurface. 

Then we considered the cosmological perturbations on superhorizon scales. We found that the effective theory on 
the brane is still the same as the case of the homogeneous and isotropic background. In particular, an anticipated 
non-trivial contribution from the spatially anisotropic part of E^^^, was found to vanish at first order in the low energy 
approximation, i.e., at 0(i?^£^). This implies that possible braneworld signatures may appear only at 0{H'^i'^) 
or higher, or second order in terms of the slow-roll parameter /H"^. Very recently, based on their low energy 
expansion method [27], and by appealing to the AdS/CFT correspondence, Kanno and Soda obtained a low energy 
effective action for the dilatonic braneworld [28]. Our result that braneworld signatures may appear only at 0{H^£'^) 
is completely consistent with their recent result [28]. 

To predict signatures specific to the braneworld scenario, we thus have to investigate the effects of 0{H^l'^). Also, 
it will be interesting to consider the other, high energy limit H'^P ^ 1 in the bulk inflaton model. Further, it is 
certainly of interest to explore the geometry in the bulk in more details, particularly the structure near the future 
Cauchy horizon of AdSs . We hope to come back to these issues in the near future. 
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APPENDIX A: PROJECTED ENERGY MOMENTUM TENSOR 

In this Appendix, we write down all the projected components of the energy-momentum tensor that appear in the 
effective 4-dimensional Einstein equations on the background metric given by Eq. (2.14). Specifically, we consider the 
quadratic part of the energy- momentum tensor, 

T^^ = da(t> db<l> - 9ah (^ff^''5c'/> da(l> + \m'(f) , (Al) 
and write down the components for the spatially homogeneous and inhomogeneous cases separately. 
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1. Spatially homogeneous background 

When is spatially homogeneous, we have 



^ab y* 1i — ^' 
n(2) a„6 _ 



b ' 



ab 



where the dot (') denotes dt and the prime (') denotes d^- 

The 4-dimensional projection of the bulk energy- momentum tensor is 



r(b,2) 



The components of this tensor are given by 



1. 



T(b,2) _ 



= 0. 



7:^'2) = ^(30'2+5<^2_3^2^2^2^5 

3.(b,2) ^ 3.(b,2)a^ ^ JL (^^'2 + ^2 _ ^2^: 



2 J,2 



(A2) 



(A3) 



(A4) 



2. Perturbation 

Here we consider the case when cj) has a spatially inhomogeneous perturbation 



IS 



The components are 



5T. 



(2) 



ah 



2d(^a4>db)H - gab{g'"^dc(t>dd5(t> + m^(t)54> 



6T. 



(2) 



ab 



qlq'l = (^(fl 6(1)' +<j)5(j)^ + m%'^(j>5(l) 



):7^(2) a & 
OJ-ab 'it 



2 



The linear perturbation of T^f^ 

(A5) 



- m b (pdtj) 5ij, 



(A6) 
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The perturbation of the 4-diniensional projection of the bulk energy-momentum tensor is 

'ST!:f qldt + {ST('J - ^5T(^))g„J . (A7) 



™(b,2) _ ^ 
d^ab - 3 



The components are 



2 

2 



^2^(b,2) ^ -(^cj)'6(j)' +^5(t>- mH'^(j)6(j)y (A8) 

APPENDIX B: SOURCE TERMS FOR E^u 

Here, wc give expUcit expressions of the terms that contribute to the source term S^v in Eq. (3.10) for iJ^i^ in the 
bulk. Again, we treat the spatially homogeneous background and the perturbation separately. 
The tensors that appear in the source term are 

p [7^(2)1 _ _^^('r(2) _ 1t(2)\„ _ 2^r(2) _ i^T(2)a„ _ln ( ^(2) _ \'r(2)\ 
^/^•^i-^ab 1 ~ 3 {)2 \ 4 J^^^'^ ^yi 3 52 « ^f^" ^■^ny-nn J'iH'^ 

^,[T^h = D^T^"/^ + 2^r£)g« , (Bl) 
where Q^va contributes to the source term only through the form i?"Qa(/i!^)- 



1. Spatially homogeneous background 

For the spatially homogeneous background discussed in Sec. Ill C, the non- vanishing components are 
P rT(2)l 1 2.,.,, 36' ,2 ,0V a<P'^ b'^^ H' 0> 



All the other components of P^i,, D"Qa{ij,v) and vanish. 

2. Perturbation 

We list the explicit expressions for the source terms under the presence of a perturbation = x{tj z)Y{x^). As 
discussed in Sec. IV, we adopt the expansion in terms of spatial harmonics Y{x^). 

The relevant coefBcients of the harmonic expansion are defined in Eqs. (4.9), (4.10) and (4.12). They are given by 

a 6' • 1 • . 1 ■ , 1 

fc /d 2 ^, 26' . 1 . 1 . , 2 ^, 



K fa -1 ^, 20' ■ i • i r , 2 „ .\ 



P. = -|^*\. (B3) 
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Q = 



la 



- — — m (px 



b' a I a 1 
0"^ a a 0'^ 

1,2 



a b' ■ , a I ■, , a 1 ■ ,, 1 a 2 , • 

a 0-^ a 26^ a 0^ 2 a 



"a2 



a 62 
0^ o 2o^ 



Qi = — 



a 



rn^ a , d b', d 1 
6 a a 3o'* a 3o- 

lav 1 



1 



1 

662' 



1 



Q2 
S = 

Si = 



fc2 



a2 L6 



+ 



k 

a . 



262' 



1 •;. a 1 

262^^" a 262 
// d 2 ■ 1 
a 36^ 

, ^ ^ ,„ m2 • a2 1 
^+a36^^^+-^'^^-^^ 

i ^ d 1 

662 Q, gj,2 

1 

362 

662"^'^' " 362 
^0' ■ , I u , 26' 4d 1 • 

-^<^x+^<^x+-^<^x--^<^x- 

'•/ ^ ^- 2 fc2 
^ 26^'^^- 36^^"^^' 



V-^^x 



362 '^'^ ~'~ QJj2 J 

d2 1 . a 1 • 6' ■ 

-^26^^^+ a 66^"^^ -66^^ ^ 
6' ■ 1 • 1 • 



>'x' 



^^-36^ 
'<^X-2 



''x- 

!''X+' 

2^- 26' 
■m (px + 



111 ■ a 1 • 1 
^^+a36^'^^ + 
1 



6^^^+ 6^' 
„ a 1 • 1 

a 36^ 



X' + 



1 



662 

fc2 1 . 1 

^66^'^^J 

3^^^_ 



--TO 0X+ -^(PX 



26^ , 2d I 



za i • 2 



(B4) 



(B5) 



APPENDIX C: GREEN FUNCTIONS FOR E^^ 



In this Appendix, we analyze the late time behavior of the retarded Green functions for Eq. (3.19) for Eu and 
Eq. (4.29) for the spatially anisotropic part E2Yij. We consider the long wavelength limit. Hence Eq. (3.19) is the 
same as that for the spatially homogeneous case, Eq. (3.22). Since the difference between these two Green functions is 
rather trivial, we mainly analyze the Green function for Ett- The Green function for E2 is given in the last subsection. 
Our analysis here is parallel to the one given in [12]. 



1. Mode functions 

We construct the Green function from a superposition of mode functions that satisfy the source-free equation. 



b{z)d, (^5.) - - 7Hdt - 12H^] ^{z, t) = 0. 



Putting ^' = Uq{z)'il)q{t), we have 



6(^)9. (^a,) + (n^ - 2)if2]u^(^) = 0, 

92 + IHdt + (n^ + 10)i/2] ^g(t) = , , 

where we have put n"^ = 0^ + 9/4 for later convenience. 

Equation (C2) can be rewritten in the standard Schrodinger from by setting Vq{z) = b~^/^{z)uq{z), 



[-—+H^V{z)\vq{z)=v:'H\q{z), 



where 



V{z) 



4 Asm\,^[H{\z\+Zo)\ 



+ H cot\i[H{\z\+ zq)\8{z) . 



(CI) 

(C2) 
(C3) 



(C4) 



(C5) 



The potential has a positive delta-function singularity at 2; = 0, i.e., it has an 'anti- volcano' form. It then follows 
that there exists no bound-state mode [17], imlike the case of the bulk scalar field itself. As the potential approaches 
9/4 at l^l 00, 71,2 > 9/4 and the spectrum is continuous. Thus, we may label the mode functions in terms of 
q = ±-^n2 — 9/4 as ^q{t, z) = Uq{z)ipq{t) with —00 < q < 00. 
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There are two independent mode functions belonging to each eigenvalue q. For a given q, we set 

1 



which satisfies the ortho-normality and completeness, 



277 



^-iqHt^-lHt 



f 



He™* 



(C6) 



(C7) 



With this choice of ipq (t) , a convenient choice for the two independent solutions is to require the outgoing- wave and 
ingoing- wave boundary conditions at |z| = oo. Specifically, 



1 m-iq) 



^ iq ^ -iq 



where Qv{C) is the Legendre function of the second kind and 2^1 [a,/?; 7, y] denotes the hypergeometric function [29], 
and ^ = cosh[iJ(|2:| + 2:0)] • Since the regions 2; > and 2: < are identical, we focus on the region 2 > in what 
follows. 



2. Construction of the Green function 

We now construct the retarded Green function introduced in Eq. (3.23) from a superposition over the modes, 

/C30 
dqGq{z;z')ijS)^S'W{t') , (CIO) 
-00 

where Gq{z,z') satisfies 



[b{z)d.(-^^d.)+{n'-2)H^]G,iz;z') 



sinh[iJ(2 + 20)] 



5{z-z'). 



(Cll) 



To construct Gq{z,z'), we need to require the Z2-symmctry, or the Neumann boimdary condition on the brane, in 
addition to the outgoing-wave boundary condition at (conformal) infinity. The Z2-symmetric mode function is given 

by 



(z) = (2) - jqu"-"^' (2) , 7g 



.(out) 



(in). 



a^uM(2) 



z=0 



Then we have 



where Wg is the Wronskian defined by 



Gg{z;z') = ^[u(°-'\z')u(g^-\z)e{z' - z) + u(°-'\z)u[^-\z')e{z - z')] , 



(C12) 



(C13) 



Wq ^ ^^"^[^^; + ^°)] (.)„(out) (,) _ ^) (.)a.n(°-) (2)] = ^ ^ . (C14) 

The behavior of the Green function at late times can be studied by analyzing poles of Gq, or equivalently zeros of 
the Wronskian, in the complex g-plane [12]. For this purpose, here, we give the 2-derivative of explicity, 

a^tt(°"*)(2) = (-i-Fig)2-5+*9^-|+*9sinh[if(2-h2o)] 



X 2F1 



''^-iq l-iq 



2 ' 2 '"^ 



,.-2 i(i-»9) 

1 — iq 



\-iq l-iq 



(015) 



19 



From this, expressions for and dzU^^^^^ follow immediately. 

Hereafter, we consider the late time behavior of the Green function in the cases HI <^ 1 and Hl'^ 1 separately. 

1. The case Hl^l 

From Eq. (C15), 

a.u(-*)(0+) « 2-H^.gsinh(g.o)(-^+ig)eo"'^^+^^ 3 ^^^t , (C16) 

and similarly for a^u('")(0+). Thus, in the limit HI < 1, 



-1/2 + iq r(l - ^g)^(| + ig) 
-1/2 -iq T{l + iq)T{l-iq) ' 



7. = tttS! ;;: , z;;;^ > (ci7) 



and the Green function Ga becomes 



G,{z-z') = {Hl)H'[inis.n\,{q^)P_^+,^{(i)P_.+,^{t) 

+ iP^^+,,m-l+^,iC')0{z' -z) + g_i+,,(e)P_ i+,,(C>(z - z'))] . (C18) 

The total Green function is given by Eq. (CIO), i.e., by integral with respect to q along the whole real axis on the 
complex g-plane. The retarded condition is guaranteed if there is no pole in the upper half of complex g-plane, which 
is in fact the case. 

For t > t', the integral (CIO) is equivalent to the contour integral over the whole region of the lower half complex 
q'-plane, or the sum of the residues of these poles. In the limit H£ <C 1, we find that Gq has poles only at 

q = -i^l±l {s = 0,1,2, ■■■). (C19) 

Hence at t > t', we may express the retarded Green function as 

oc 

G{t,z;t',z') = {Hl)9{t-t')J2Pn{OPnie)a-^^''+^Ht)a'^+Ht'), (C20) 

n=0 

where P„ (z) is Legendre polynomial of order n. In particular, after a sufficient lapse of time, we may approximate the 
Green function the n = term, i.e., by taking account of only the pole with the smallest imaginary part q = —i/2. 
This gives 

G{t, z; t', z') « (HI) 9{t - t') a-'^{t)a'^{t') . (C21) 

Thus, G oc a~'^ at sufficiently late times. This gives rise to the radiation-like behavior of Eu, in accordance with our 
expectation. 

2. The case i?Z » 1 

In this case, we have 

7<, = ^^^2^-eo^-, (C22) 

and the structure of Gq{z,z') is more complicated than the case Hi <^ 1. Nevertheless, it is possible to locate the 
poles of Gq. We find the poles at 

9 = -^, -^s (s = l,2,3,---). (C23) 

Thus at sufficiently late times, the Green function behaves as 

G{t, z; t', z') « e{t - t') a-^{t)a^{t') . (C24) 

Apart from the absence of the factor HI, this is the same as that in the case H(. <C 1, Eq. (C21). Hence the 
radiation-like behavior is realized again. Of course, this is also expected, because the transverse-traceless nature of 
the source-free E^i, is independent of the energy scale. 
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3. Green function for E2 

The Green function G2{t,z;t',z') for E2 satisfies 

'b{z)d, (^3^) - dt - ^Hdt - 2ir2] G2{t, z; t', z') = -6{t - t')5{z - z')b{z), 



(C25) 



in the long wavelength limit. In this case, only the difference from the Green function for Eft appears in the temporal 
part of the mode functions. For an eigenvalue q, with q'^ = — 9/4 as before, the ^;-component of the equation is 
exactly the same as Eq. (C2), while the f-component becomes 



= 0. 



The solution is related to tjjqit) for Eu as 'ip2,g{t) = {t)'tpq{t) . Hence, it follows that 

G2it,z;f',z') = a'^it)G{t,z;t'z')a-~^it') , 



(C26) 



(C27) 



and we can use all the results in the previous subsections. In particular, in the limit H£ <C 1, the Green function at 
t > t' is expressed as 



G2{t,z;t',z') = {Hl)9{t-t')J2Pn{OPn{aa-^^^'Ht)a^^^{t'). 



(C28) 



n=0 



Also, at sufficiently late times, we have G2 (x a ^, irrespective to HI. 
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